Search for low lying dipole strength in the neutron rich nucleus 26 Ne 
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Coulomb excitation of the exotic neutron-rich nucleus 26 Ne on a na *Pb target was mea- 
sured at 58 A.MeV in order to search for low- lying El strength above the neutron emission 
threshold. Data were also taken on an na *Al target to estimate the nuclear contribution. 
The radioactive beam was produced by fragmentation of a 95 A.MeV 40 Ar beam deliv- 
ered by the RIKEN Research Facility. The set-up included a Nal gamma-ray array, a 
charged fragment hodoscope and a neutron wall. Using the invariant mass method in 
the 25 Ne+n channel, we observe a sizable amount of El strength between 6 and 10 MeV. 
The reconstructed 26 Ne angular distribution confirms its El nature. A reduced dipole 
transition probability of B(Ei) = 0.49 ± 0.16 e 2 fm is deduced. For the first time, the 
decay pattern of low-lying strength in a neutron-rich nucleus is obtained. The results are 
discussed in terms of a pygmy resonance centered around 9 MeV. 



1. Introduction 

Giant Resonances are a general feature of nuclei and their properties give us a handle 
on the macroscopic and microscopic behavior of nuclear matter. These modes have been 
extensively studied in stable nuclei over the last 50 years and the recent inception of 
Radioactive Ion Beam facilities opens the opportunity to extend these investigations to 
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exotic nuclei. Far from stability new modes are expected to appear. In particular, the 
dipole response of neutron-rich nuclei may exhibit strength at energies lower than the 
standard Giant Dipole Resonance (GDR) often depicted as the oscillation of a deeply 
bound core against a neutron skin giving rise to a so-called pygmy resonance. Recent 
studies on Oxygen [lj and Tin [2] isotopes have given the first experimental indications 
for such an effect, while at the same time calling into question its nature. 

The theoretical approach is often based on mean field calculations. Relativistic random 
phase approximation (RRPA) and quasi-particle RRPA (QRRPA) calculations have been 
carried out by Cao and Ma in 26 Ne [3] which predict that below 10 MeV excitation energy, 
almost 5% of the Thomas-Reiche-Kuhn (TRK) energy weighted sum rule is exhausted by 
strength centered around 8 MeV. This region of energy is located between the one neutron 
and the two neutron emission thresholds. Thus, in order to investigate this prediction, we 
performed Coulomb excitation of 26 Ne at intermediate energies on a lead target and used 
the invariant mass method to reconstruct the B(E1) strength from the 26 Ne— » 25 Ne* + n 
reaction. 

2. Experimental details 

The experiment was performed at the RIKEN Accelerator Research Facility. A sec- 
ondary 26 Ne beam was produced through fragmentation of a 95 A. MeV 40 Ar primary 
beam on a 2-mm- thick 9 Be target. The 26 Ne was separated by the RIKEN Projectile 
Fragment Separator (RIPS) [4J. Particle identification was unambiguously performed by 
means of the time-of-flight (TOF) and the purity was 80%. The 26 Ne beam of intensity 
~ 5.10 3 pps and incident energy 58 A. MeV, was tracked with two parallel-plate avalanche 
counters providing incident angle and hit position onto the target. It then impinged 
alternatively on a 230 mg/cm 2 na *Pb and a 130 mg/cm 2 na *Al target. Data obtained 
with aluminum target are used in the following to estimate the contribution of nuclear 
excitation to the data. 

The outgoing charged fragments were measured using a set of telescopes placed at 
1.2 rn downstream of the target. They consisted of two layers (X and Y) of 500 /xm 
silicon strip detectors (SSD) with 5 mm strips which provided an total energy resolution 
of 1.5 MeV (FWHM). The last layer used 3-mm-thick Si(Li) from the charged-particle 
detector MUST [5], provided a 9 MeV (FWHM) resolution on the remaining energy (E). 
Unambiguous mass and charge identification of all projectile like fragments was obtained 
using the E-AE method. 

In-beam gamma rays were detected using a 47r-gamma-detector, DALI2 [6] , which con- 
sists of 152 Nal detectors placed around the target. For 1.3 MeV gamma-rays, its measured 
efficiency is approximately 15% with an energy resolution of 7% (FWHM). The Doppler 
corrected gamma energy distribution obtained in coincidence with the 25 Ne isotope al- 
lows us to identify the gamma decay from the adopted 1702.7(7) keV and 3316.4(11) keV 
excited states. In addition we observed a 2063 ± 67 keV gamma-ray which we related to 
the 2030 ± 50 keV excited state, only seen up to now by transfer reaction [7j. 

The hodoscope for neutron detection was an array of 4 layers of 29 plastic rods each, 
placed 3.5 m downstream of the target. Each layer was composed of 13 [2.1 mx6x6 cm 2 ] 
and 16 [1.1 mx6x6 cm 2 ] rods, arranged in a shape of a cross. Its total intrinsic efficiency 
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for the detection of one 60 MeV neutron was calculated to be ~ 25%. Finally 29 thin 
plastic rods covered the front face of the wall in order to veto charged particles as well as 
to provide an active beam stopper. The neutron position is determined with an error of 
±3 cm and the energy, from TOF information, with a 2.5 MeV (FWHM) resolution for 
the neutrons of interest. 

3. Results 

We have performed a simulation of the experimental setup using the Geant 3 pack- 
age [8]. In order to check the reliability of the analysis, we built the elastic scattering 
angular distribution of 26 Ne on na *Pb at 55 A. MeV. It is in good agreement with a DWBA 
calculation based on a 20 Ne+ 208 Pb at 40 A. MeV optical potential from [9]. We empiri- 
cally generated optical potential parameters for the 26 Ne+ 27 Al reaction and tested them 
by comparing with the experimental elastic scattering. 

Using the invariant mass method, the excitation energy of an unbound state in the 
nucleus A X decaying to a state in A_1 X can be expressed by: E* = E Ie \ + S n + ^-E 7i , 
where E re \ is the relative energy between neutron and A_1 X and S n the one neutron 
emission threshold. The gamma detection efficiency was not high enough to apply an 
event-by event reconstruction technique. To extract the excitation energy spectrum, we 
hence used a method based on adequate subtraction of relative energy spectra built for 
these two kind of event. 

The method can be illustrated in the schematic case where the daughter nucleus A-1 X 
has only one excited state below its neutron threshold. The excitation energy spectrum 
can be decomposed in the sum of two contributions: the decay of A X to the ground-state 
(gs) of A_1 X and the decay of A X to the excited state (I) of A_1 X. 

The first contribution can be obtained by subtracting from the inclusive relative energy 
spectrum [£ , re i( A_1 X, n)] the relative spectrum of gamma-coincidence events [£ , re i( A_1 X,n)| 7 ] 
divided by the gamma detection efficiency e, and shifted " c — >" by the one neutron emission 
threshold S n i.e.: [E*] 9S = [E Iel ( A -%n)] - [E Iel ( A '% n)\ 7 ] /e - S n . 

The second contribution is simply the relative energy spectrum of gamma-coincidence 
events shifted by the gamma energy: [-E 1 *!^ = |_-^ rel ( A_1 -^' ra )l7 1 / e $n "—>■ E 7 . Fi- 
nally, the excitation energy spectrum for the 25 Ne+n decay channel is given by: \_E*~\ = 
[E*] 9S + [E*] m . The method was tested by simulation in this schematic case, and also 
in the realistic case taking into account the experimental decay scheme of 25 Ne. 

The excitation energy spectra reconstructed for the 25 Ne+n decay channel obtained with 
the nat Pb and nat Al targets are represented on Fig. [TJ Note that above S2 n = 9.8 MeV, the 
decay of 26 Ne is expected to occur mainly by 2 neutron emission. Between 8 and 10 MeV, 
a sizable amount of cross-section is observed for both targets. In intermediate energy 
inelastic scattering with a heavy target such as raa *Pb, the Coulomb dominance of the 
El excitation is well-known. The contribution of possible E2 excitation to the spectrum 
obtained with the lead target has — in a first step — been determined using data taken 
with the aluminum target and the coupled channels ECIS 97 [10J code. Assuming simple 
collective vibration mode with equal nuclear and Coulomb deformation lengths, the E2 
deformation parameters were extracted from the measured cross-section with the na *Al 
target (oai = 9.1 ± 2.3 mb). The L = 2 cross section in lead was then calculated using 
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Figure 1. Left: Excitation energy distribution in 26 Ne in the nat Pb target. The shaded 
area is a tentative Lorentzian fit. Right: Same as previous but for the nat Al target. 



the deformation lengths extracted in the previous step. We obtained (Tp^ 2 = 17.9 ± 
4.3 mb. After subtraction of the L = 2 contribution, the resulting (Jp^ 1 = 48.5 ± 4.8 mb 
cross section corresponds, using ECIS 97, to a Coulomb deformation parameter of fie — 
0.087 ± 0.05 which led to a B(El) = 0.55 ± 0.05 e 2 fm 2 via the relation B(E1; 0+ -» 
1~) = {~§^Z p e RcPc =1 ) 2 with R c the Coulomb radius. This value of reduced transition 
probability correspond to 5.5 ± 0.6% of the Thomas-Reiche-Kuhn energy weighted sum 
rule for an excitation energy of 9 MeV. We estimate the error due to the choice of optical 
potential by performing the same analysis using now parameters from the 40 Ar+ 208 Pb 
reaction at 40 A. MeV [IT]. We obtained B(El) = 0.60 ± 0.06 e 2 fm 4 i.e. 6.0 ± 0.6% of 
the TRK, in good agreement with the values previously extracted, which demonstrates 
that we are not strongly sensitive to the choice of the optical potential for the reaction 
on lead. 

Due to the high granularity and the good resolution of the present setup, it is possible 
to reconstruct the scattering angular distribution for 26 Ne on the ™ a *Pb target which is 
represented in Fig. [21 Hence, we have used a second method to extract the El excitation 
which relies on a multipole decomposition analysis of this angular distribution. The 
L = 1 and L > 2 angular distributions (dashed and dotted lines) were obtained from 
simulation based on ECIS 97 angular distribution calculated for E* = 9 MeV. The data 
were fitted with a linear combination of the two distributions. The results of the fit give us 
B(E1) = 0.49±0.16 which corresponds to 4.9±1.6% of the TRK, again for an excitation 
energy of 9 MeV. If we suppose now that the remaining part of the contribution is due to 
L = 2 excitation we can extracted a B(E2 f ) = 49 ± 8 e 2 fm 4 for this structure. 

The two results from the two different methods are in agreement but since the multipole 
decomposition of the angular distribution relies only on the data from lead target, the 
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Figure 2. Result (solid line) of the mul- 
tipole decomposition of the experimen- 
tal angular distribution (circle) of the 
26 Ne angular scattering onto lead tar- 
get, by L = 1 (dashed line) and L = 2- 
type (dotted line) distributions. 



second value of B(E1) = 0.49 ± 0.16 e 2 fm 2 will be retained. The results obtained from 
our experiment concerning the El transition are compared in the following to theoretical 
calculations. 



4. Comparison with theory 

In the introduction we already presented results from Cao and Ma [3]. Using the 
relativistic QRPA framework and the response function formalism they predicted a pygmy 
resonance centered around 8.4 MeV and which exhausts 4.5% of the TRK, which is close 
to our experimental values. Another calculation has been performed by Khan et al. |12j . 
It is based on effective SGII Skyrme interactions and is performed in the spherical QRPA 
framework. It predicts a redistribution of the strength at low energy centered around 
E* = 11.7 MeV exhausting ~ 5% of the TRK. Two other preliminary calculations has 
been performed in the deformed QRPA framework using Gogny forces [13] and in the 
deformed relativistic QRPA framework [TJ]. Both predict a redistribution of the strength 
at low energy, centered around 10.7 MeV and 7.5 MeV respectively. The first calculation 
also predicts that only ~ 1% of the TRK should be exhausted. Preliminary shell-model 
calculations have also been performed by Nowackiet al. [15] who predict a state at 9.3 MeV 
exhausting ~ 5% of the TRK. All these theories agree on the presence of a 1~ structure 
a low excitation energy, compatible with our experimental result but they disagree on its 
nature: the amount of strength differs as well as as well as its collective or single particle 
nature. 



5. Decay of pygmy resonances in 26 Ne 

Our reconstruction method for the excitation energy allows us to extract for the first 
time data on the decay of pygmy resonance of neutron-rich nuclei. We present the exper- 
imental branching ratios towards the various states of 25 Ne in Tab. [TJ The clear difference 
between the branching ratios obtained with lead and aluminum targets proves that states 
of different nature have been excited. For comparison, we also performed a statistical 



6 



Table 1 

Experimental branching ratio compared to statistical model for a given 26 Ne transition 
multipolarity to a given excited state in 25 Ne compared to experiment. 



Final 


25 Ne state 


Experiment 




Statistical decay 




Label 


J n 


Pb 


Al 


L = 


1 L = 2 


L = 3 


(gs) 


1/2+ 




< 10% 


40% 


28% 


22% 


(I) 


5/2+ + 3/2+ 


66% ± 15% 


95±\ 5 % 


55% 


67% 


75% 


(II) 


3/2" 


35% ± 9% 


5±|% 


5% 


4% 


3% 



decay calculation for L — 1, 2, 3 states using the CASCADE code [16], assuming spins and 
parities of 25 Ne states as listed in Tab. [TJ It clearly shows that the decay is not statistical, 
which is not surprising for a light nucleus. No predictions of the direct decay of pygmy 
states yet exist from the previously mentioned microscopic models. Future comparisons 
with our data should be a strong test for these models. 

6. Conclusion 

We performed the Coulomb excitation of 26 Ne in order to measure its low lying dipole 
strength below 10 MeV excitation energy, using the invariant mass method. We extract 
an El strength value of B(E1) = 0.49 ±0.16 e 2 fm 2 between the one- and the two-neutron 
emission threshold as well as the corresponding decay pattern. Our results are compatible 
with various theoretical predictions. Future comparison of the extracted decay pattern 
with theoretical models may allow us to elucidate whether the strength is of single particle 
nature or rather due to a collective pygmy resonance. 
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